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ABSTRACT 

We consider essential mechanisms for orbit-shrinkage or "hardening" of compact bina- 
ries in globular clusters to the point of Roche-lobe contact and X-ray emission phase, 
focussing on the process of collisional hardening due to encounters between binaries 
and single stars in the cluster core. The interplay between this kind of hardening and 
that due to emission of gravitational radiation produces a characteristic scaling of the 
orbit-shrinkage time with the single-star binary encounter rate 7 in the cluster which 
we introduce, clarify, and explore. We investigate possible effects of this scaling on 
populations of X-ray binaries in globular clusters within the framework of a simple 
"toy" scheme for describing the evolution of pre-X-ray binaries in globular clusters. We 
find the expected qualitative trends sufficiently supported by data on X-ray binaries 
in galactic globular clusters to encourage us toward a more quantitative study. 

Key words: globular clusters: general - scattering - stellar dynamics - binaries: 
close - stars: low-mass ~ stars: horizontal branch - X-rays: binaries 



1 INTRODUCTION. 



It is well-known that globular clusters contain far more than 
their fair share of compact X-ray binarie s per unit stel- 
lar mass, compared to their h ost galaxies iVerbunt fc Hud 
Il987l. IVerbunt fc LewTnl |2004 . The enhancement factor is 
100 in the Milky Way and M31 ifVerbunt fc Lewin 20*01 
iPooley et.aD l2003!l. and possibly much higher in elliptical 
:alaxies, as recent Chandra observations have suggested 



a low-mass normal companion, although such binaries can 
be low-luminosity X-ray sources. However, general consider- 
ations on the dynamical formation of close binaries do apply 
to these as well; indeed, the latter binaries ar e now widely ac- 
cepted as e volutionary products of LMXBs ivaa den Heuvel| 
Il99llll992l). 



lAngelini et.ai]l200lL iPoolev et.af]l2003^ . The origin of this 
overabundance of close binaries has been realized for some 
thirty years now to be the dynamical formation of such bina- 
ries — through tidal capture and/or exchange interactions 
— which can proceed at a very significant rate in dense 
cores of globular clusters (hencef orth GCs) because of the 
high stellar-encounter rates there (lHutll985l . lHut fc Verbun3 



ll983l.lHut et!al]|l992l) . but whose rate is negligible over the 
rest of the galaxy, where the stellar density is low by com- 
parison. The GC X-ray binaries that we shall be mainly 
concerned with in this work are those which are powered by 
accretion onto compact stars. These can be either (a) low- 
mass X-ray binaries (henceforth LMXBs), containing neu- 
tron stars accreting from low-mass companions, or, (b) cata- 
clysmic variables (henceforth CVs) , containing white dwarfs 
accreting from low-mass companions. Accordingly, we shall 
not explicitly consider here binaries which contain either 
(a) two "normal" solar-mass stars, one or both of which are 
coronally active, or, (b) recycled neutron stars operating as 
rotation-powered millisecond pulsars, with a white-dwarf or 



In tidal-capture formation of a compact-star binary, the 
compact star (neutron star/white dwarf) passing close to a 
normal star dissipates its kinetic energy significantly by cre- 
ating tidal deformation in the latter star, and so becomes 
bound to it. In the exchange process of formation, a com- 
pact star replaces one of the stars of an existing binary sys- 
tem of two normal stars during a dynamic al encounter iHud 
Il985l iHut fc Verbundligil ISpitzeJll987l) . After formation 
in such an encounter, the compact-star binary continues to 
undergo stellar encounters in the dense cores of GCs, and 
it is on one particular effect of the continuing encounters 
that we focus in this paper. In the mid-1970s, it was real- 
ized that a major effect of the binary-single star encounters 
would be to extract energy from a given binary, making it 
more tightly bound or harder, and giving this energy to the 
mot ion of the single stars in the GC, thus "heat ing" the clus- 
ter jHeggijll975l . lSpitzedll987l.lHut et.al.lll992l) . We can call 
this effect collisional hardening of the compact-star binary, 
which makes the binary's orbit shrink at a rate higher than 
that which would obtain if it were not subject to the above 
stellar encounters, i.e., if it were not in a GC. The latter rate 
is believed to be determined by a combination of two pro- 
cesses, VIZ., (a) emission of gravitational radiation and (b) 
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magnetic braking. In the former process, a compact binary 
emits gravitational radiation and so loses e nergy and angula r 
momentum, which makes its orbit shrink l)d'Invernolll992il . 
The latter mechanism is envisaged as follows. The low-mass 
companion to the compact star has a significant magnetic 
field, and also has its rotation tidally coupled or "locked" 
to that of orbital revolution. The companion drives a wind, 
which carries away angular momentum at an enhanced rate 
because the magnetic field enforces co-rotation of the wind 
out to a radius considerably larger than that of the star, and 
this angular momentum ultimately comes from the orbit be- 
cause of the above tidal lo cking, thus making the orbit shrink 
JVerbunt fc Zwaan|[l98lh . We discuss in this paper the rel- 
ative roles of the above mechanisms for binary hardening, 
particularly the role of collisional hardening vts-d-vis that 
due to gravitational radiation, indicate and clarify a scal- 
ing that naturally emerges from this interplay, and briefly 
suggest possible observational signatures of this scaling. 

In Sec. 2, we discuss the hardening of compact bina- 
ries by the three mechanisms discussed above, bringing out 
the particular role of collisional hardening. We show that 
a characteristic scaling of the orbit-shrinkage time of such 
binaries with an essential GC parameter emerges because of 
the interplay between collisional hardening and that due to 
gravitational radiation. In Sec. 3, we explore possible obser- 
vational signatures of this scaling. We sketch a very simple 
"toy" scheme for describing the evolution of compact X-ray 
binaries in GCs, and indicate a possible signature of the 
above scaling within the bounds of this toy model. We show 
that current data on X-ray binaries in GCs are consistent 
with this signature. We discuss our results in Sec. 4, explor- 
ing possible lines of future inquiry. 



2 HARDENING OF COMPACT BINARIES. 

We consider the shrinking of the orbital radius of a compact 
binary by the three mechanisms introduced above. Consider 
gravitational radiation flrst. The rate at which the radius a 
of a binary decreases due to this process is given by (see, 
e.g., d'Inverno 1992 and references therein): 

-3 



fflGw = aawci 



aaw ~ -12.2Mmxmc Re/Gyi (1) 



In this equation, mx and tUc are respectively the masses of 
the compact star and its companion in solar masses, M = 
mx + rric is the total mass in the same units, and the orbital 
radius a is expressed in units of solar radius. We shall use 
these units throughout this paper. 

Now consider magnetic braking. The orbit shrinkage 
rate due to this pr ocess is given in the ori ginal Verbunt- 
Zwaan prescription iVerbunt fc Zwaanlll98l[) as: 

mx 

where Rc is the radius of the companion. We give below 
further discussion on this mechanism. 

Finally consider the rate of orbit shrinkage due to col- 
lisional hardening, which is given by (IShulllll97gD . 

3 

dc = ac7a^ ac ~ -2.36 x 10~^ -Ro/Gyr (3) 

mxmc 

Here, mac is the mass of the single stars in the GC which 
are undergoing encounters with the binary. In this introduc- 
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Figure 1. Relative orbit shrinkage rates -a/a due to gravitational 
radiation, magnetic braking and collisional hardening, shown as 
functions of the binary separation a. Also shown is the total 
shrinkage rate. Value of 7 as indicated. 



tory work, we assume mac to be a constant, representing 
a suitable average value for a GC core, which we take to be 
mac ~ O.6M0. 

The parameter 7 is a measure of the encounter rate 
between a given binary and the background of single stars 
in the core of the GC: it is a crucial property of the GC for 
our purposes, so that we shall use it constantly in this work. 
It scales as 7 oc pc/vc with the (average) core density pc of 
the GC, and the velocity dispersion Vc of the stars in the 
core. Follo wing the convention ofte n employed in the GC 
literature JVerbundl200l lHudll985l) . we can, in fact, define 
this parameter as: 



7 ■ 



(4) 



Then the unit of 7 is « 6.96 x IO^Mq^Rq'* sec, corresponding 
to the units of pc and Vc commonly used in the GC literature, 
namely, Mqpc"'' and km sec~^ respectively. In these units, 
values of 7 generally run in the range ~ 10^-10'' (see below). 

In this work, we take the mass of the compact star to 
be mx ~ 1.4Mq and that of the companion mc — 0.6Mq, 
the latter corresponding to a typical average mass of stars in 
a GC (see above) . According to the mass- radius relation for 
low-mass main-sequence stars, the radius of a main-sequence 
companion will then be Rc ~ O.GRq. Furthermore, we con- 
sider only circular orbits in this work, returning to this point 
in Sec. 4. 

The total rate of orbit shrinkage due to the combination 
of the above mechanisms is given by: 



a = dew + fflMB + etc 



(5) 



We emphasize that the first two terms in Eq. @ are always 
operational, irrespective of whether the binary is in a GC 
or not, and it is the relative effect of the third term, which 
represents the effects of the encounters in a GC core, that 
we wish to study here. The interplay between the first and 
the t h ird te rm was investigated in a pioneering study by 
IShulJ (Il97fll). before the magnetic b raking mechanism was 
postulated iVerbunt fc ZwaanlUgSlTl . 
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Note first that the three terms have different regions 
of dominance, as shown in Fig. Coiiisionai hardening 
dominates at iarge vaiues of the orbital separation a, i.e., 
for wide binaries, while hardening by gravitational radia- 
tion and magnetic braking dominates at small a, i.e., for 
narrow binaries. Between the latter two, magnetic braking 
dominates at the smallest orbital separations, if we adopt 
the original Verbunt-Zwaan (henceforth VZ) scaling for it 
(see below). The relative orbit shrinkage rate a /a thus scales 
as a at large orbit separations, passes through a minimum 
at a critical separation ac where the gravitational radiation 
shrinkage rate, scaling as a /a ~ takes over from colli- 
sional hardening, and finally rises at very small separations 
as a/a ~ due to VZ magnetic braking. The change- 
over from gravitational radiation shrinkage to that due to 
magnetic braking occurs at a radius am « ac. These two 
critical radii are easily obtained from Eqs. Q, ©, and 
and are given by 



1/5 

= "gw"c 7 



1/5 -1/5 



aiviB 

CfGW 



(6) 



Note the scaling Oc oc 7"^''^, which is crucial for much of 
our discussion here, as we shall see below. The critical or- 
bital separation Oc varies in the range ~ (5 — 12)7i0 for 
the canonical r ange of values of the above GC parameter 
JVerbuntlbOO^) 7 ~ 10^ - 10^ in the above units JShulJ 
119791) . 

The relevance of this to close compact-star binaries in 
GCs is as follows. When such a binary is formed, its orbital 
separation in most cases is such that the low-mass com- 
panion is not in Roche lobe contact, since the Roche-lobe 
radius has to be Rj^ ~ O.6-R0 or less for this to happen 
for a typical low-mass main sequence or subgiant compan- 
ion of mass ~ O.6M0 (see above). Mass transfer does not 
occur under such circumstances, so that such binaries are 
pre-LMXBs or pre-CVs, and we can call them by the gen- 
eral name pre-X-ray binaries, or PXBs for short. It is the 
above orbit-shrinkage or hardening process that brings the 
companion into Roche-lobe contact, so that mass transfer 
begins, and the PXB turns on as an X-ray binary (LMXB 
or CV), or XB for short. 

Depending on the initial separation ai of the binary, 
some or all of the above processes can thus play signifi- 
cant roles in shrinking it to the point where mass trans- 
fer begins. Recent numerical simulations suggest that tidal- 
capture binaries are born with orbital radii (or semi-major 
axes) in the range 1 < at/RQ < 15 for main-sequence 
(henceforth ms) or early subgiant companions, and in the 
range 40 < aj/Rp, < 100 for h orizontal-branch compan- 
ions jPortegies Zwart et.al.lll997ifl . For binaries formed by 
exchange encounters, the orbital radii are generally expected 
to be somewhat larger than those for corresponding tidal bi- 
naries with identical members. Thus, both collisional hard- 
ening and gravitational radiation are expected to play major 
roles in the orbital shrinkage to Roche-lobe contact for most 
of the PXBs in GCs, whether dynamically formed or pri- 
mordial. 

Although we have included magnetic braking as above 
for completeness, its role in hardening of PXBs into XBs ap- 
pears to be rather insignificant, at least for the VZ scaling 
adopted above. This is evident from the fact that there is lit- 
tle change in any of the results described here whether mag- 



netic braking is included or not. For the VZ scaling, this is 
easy to understand. With a steep increase at small a, this ef- 
fect is significant only at very small orbital separation, when 
the PXB has already come into Roche-lobe contact and be- 
come an XB. Thus, this process may well be significant in 
the fu rther orbital evolution of the XB a s mass-transfer pro- 
ceeds ivan den Heuvel|ll99ll Ivan den H cuvcl 19£^, but not 
in the PXB-hardening process under study here. 

Actually, further study of magnetic braking since the 
original VZ formulation has revealed many interesting 
points. The nature and strength of this effect very likely 
depends on the mass and evolutionary state of the compan- 
ion star. For example, magnetic braking may become totally 
ineffective for very low-mass companions with nic ~ O.SMq 
or less, as these stars are fully convective. Whereas a signif- 
icant convective envelope is necessary for strong magnetic 
braking, "anchoring" of the magnetic field in a radiative 
core is also believed to be essential for it, and it is argued 
that the effect w ould basically vanish when the star becomes 
fully convective dSoruit fc Ritteilll983l iPodsiadlowski et.alJ 
l2002h . Indeed, this forms the basis for the standard expla- 
nation for the period gap in CVs (van den Hcuvcl 1992 and 
references therein). Further, studies of the rotation periods 
of stars in open clusters have suggested that magnetic brak- 
ing may be less effective than that given by the VZ prescrip- 
tion: this has been modelled in recent literature by either (a) 
the VZ scaling as above, but a smaller numerical constant 
than that given above, or, (b) a "saturation" effect below 
a critical value of a, wherein the scaling changes from the 
VZ ~ a"** scaling of Eq. ^ to a much slower ^ sca l- 
ing below this critical a- value (Ivan der Sluvs et.alj I2OO5II . 
While these modifications are of relevance to XB evolution, 
it does not appear that they can alter the PXB-hardening 
results described here in any significant way. Accordingly, 
we shall not discuss magnetic braking any further here, and 
keep this term in the complete equations only to remind 
ourselves that it is operational, in principle, for companions 
with rric ^ 0.3Mq. 



2.1 An Interesting Scaling 

Interplay between collisional hardening and gravitational- 
radiation hardening near the above critical orbital separa- 
tion ac produces a characteristic scaling, which we now de- 
scribe. Consider the shrinkage time tpxb of a PXB from 
an initial orbital separation ai to the final separation a/ 
corresponding to Roche-lobe contact and the onset of mass 
transfer, given by: 



TpxB{ai,^) 



da 



daw + duB + dc 



(7) 



For given values of stellar masses, tpxb scales with the GC 
parameter 7 introduced above as 

rpxfl-7"'^'- (8) 

The scaling is almost exact at high values of 7, i.e., 7 > 10*, 
say, there being a slight fall-off from this scaling at low 7's. 

How does this scaling arise? To see this, consider first 
the qualitative features of the integrand on the right-hand 
side of Eq. ||7J, i.e., the reciprocal of the total shrinkage rate 
at an orbital separation a, which we denote by C{a), and 
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Figure 2. Integrand C{a) in Eq. shown as function of orbital 
separation a, with values of ■y as indicated. 

which is displayed in Fig. |5| It is sharply peaked at a ~ Oc: 
indeed, the peak would be exactly at the above critical sep- 
aration Oc but for the effects of magnetic braking, as can be 
readily verified. Since the latter effects are not important in 
the range of a-values relevant for this problem, as explained 
above, we can get a good estimate of the actual result by 
considering only gravitational radiation and collisional hard- 
ening. Because of this dominant, sharp peak in (^{a), most of 
the contribution to the integral, i.e., to tpxb, comes from 
there, provided that the integration limits (a;, a/) are such 
that all or most of the peak is included. We assume for the 
moment that this is so, and return to a discussion in the 
next subsection of what happens when this condition fails. 

Under the above circumstances, we can immediately 
give a rough estimate of tpxb, which is the area under 
the curve in Fig. |5| as rpx b ~ 2ac x (maximum value 
of the above integrand). This maximum value is simply 
l/{2aGwa~'^) if we neglect magnetic braking, since the 
gravitational radiation term equals the collisional harden- 
ing term there, as explained above. This gives tpxb ~ 
a^^a\, which, with the aid of Eq. JHJl, yields tpxb ^ 

'-^Gw^^c*''^')'"*^^- This is the basic reason for the scaling 
given by Eq. Q. 

An exact evaluation of the integral in Eq. Q , with the 
magnetic braking term neglected, confirms this, as expected 
and as detailed in Appendix A. The exact result is: 

TPXB = a-:^^a-*/%-^/^[/(6,) - J(60]. (9) 

Here, 6 is a dimensionless orbital separation defined by 
h = a/ttc, and the integral I{x) is given in Appendix A. 
As I{x) has only a logarithmic dependence on x under these 
circumstances, the basic scaling is tpxb ~ 7"''^^, as above. 
It is this basic scaling that leads to the essential behavior of 
the shrinkage time tpxb discussed in this paper. 

2.2 Breakdown of Scaling? 

When would the above scaling break down, and why? A sim- 
ple answer is clear from Fig. |21 this would happen when the 
integration limits {ai,af) are such that all or most of the 
above peak in (^(a) is not included. For the present problem, 



this basically reduces to an upper bound on a/, since Oi is 
normally large enough to ensure that the region of integra- 
tion in Fig. 121 extends well into considerably larger values of 
a beyond the peak. When a/ becomes so large as to exceed 
Qc, the region of integration is severely curtailed from the left 
in Fig. 121 so that most of the peak's contribution is missed, 
and the above scaling breaks down. We might think that 
such a situation would arise when the low-mass companion 
in the PXB is an evolved, horizontal-branch star, which has 
a much larger radius than a ms/subgiant companion of the 
same mass, and so would be expected to come into Roche- 
lobe contact at a much larger value of Rl, say 5 — WRq, 
and correspondingly larger values of a/. But such binaries 
are not relevant to our discussion here, since the lifetimes 
(~ lO'^ y) of such horizontal-branch stars are too short to be 
of significance to the long binary-hardening timescales under 
consideration here. Thus, this possibility is not of practical 
importance here. 

However, there is a situation in which this scaling is not 
relevant, not because it breaks down, but, rather because we 
move into a region of 7- values where tpxb computed in the 
above way exceeds the expected main-sequence lifetime Tc 
of the low-mass ms/early-subgiant companion. Under these 
circumstances, the companion starts evolving into a giant 
and rapidly fills its Roche lobe, for essentially any value that 
a is likely to have at that stage. This is formally equivalent 
to saying that Tpx b saturates at a value Tc in this range of 
7. We return to this point below. 

2.3 Shrinkage Time 

We now calculate the exact variation of the shrinkage time 
TpxB with the encounter-rate parameter 7 introduced ear- 
lier, keeping all terms in Eq. ||7J. For this, we need to specify 
the initial and final values, Oi and a/, of the orbital sepa- 
ration. We adopt a/ ~ 1.94_Rq for ms/subgiant compan- 
ions corresponding to Roche-lobe contact, when the radius 
of the Roche-lobe Rl of the companion becomes equal to 
the radius of the companion itself. This translates into the 
above value of the orbital separation a/ by the well-known 
Paczyiiski relation: 

(m \ ^/^ 

corresponding to Rl being equal to companion radius « 
0.6i?Q for a companion mass ~ O.GMq. 

In general, a; will have a value which is within a possi- 
ble range (ar™,ar''^), which is indicated in Tabled This 
range depends on the formation-mode of the binary, and 
also on the evolutionary status of the companion. The for- 
mer has two possibilities, namely, (a) the binary is primor- 
dial. I.e., it was already a binary when the globular cluster 
formed, or, (b) it formed by tidal capture or exchange inter- 
actions in the dense core of the globular cluster. The latter 
has also two basic possibilities, namely that the companion 
is either (a) a ms/subgiant, or, (b) a horizontal-branch star, 
as explained earlier. As explained above, however, the short 
lifetimes of ho rizontal-branch |Kipponlian & Woigcrt 1993, 
IClavtonllTgeil stars compared to the timescales of harden- 
ing processes under study here make it clear that they are of 
little importance in this problem, and we shall not consider 
them any further in this work. The ranges of Ui adopted in 
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Table 1. Distribution functions f{ai) and range of initial orbital 
separations ai of compact-star binaries in globular clusters 

Type of compact Range of initial Form of distribution 

binary radius Oi function ficii) 



10^ 



Dynamically formed 
compact star binary 
with ms or subgiant 
companion 



2Rq - 5OR0 



f{ai) = constant, 
f{ai) ~ ai, 
Gaussian in 
with 11 = G.ORq and 
<T = 12.7R0 



Primordial compact 
binaries 



2Rq - 500i?Q f{ai) . 



various cases are detailed in Table [H and ar e taken from 
current literature jPorteeies Zwart et.al.ll^997^ . 

Since a GC has a distribution of UiS, we wish to study 
how TpxB{'y,a.i) averaged over such a distribution scales 
with 7, since both of these represent overall properties of the 
cluster. To this end, we define a suitable average shrinkage 
time as: 



r(7) = (tpxb) 



TpxB{'y,ai)f{ai)dai 



(11) 



where /(fli) is the normalized distribution of a; in the range 
(a7"", a7*"^). For this distribution, some indications and 
constraints are available, as follows. For primordial binaries, 
the distribution f{a) oa 1/a (corresponding to a flat cumula- 
tive distribution in In a) is well-established iKraicheva et.alJ 



Il978r) . For tidal capture/exchange binaries, results from 
numer ical simulations like those of IPorteeies Zwart et.al.l 
il997l) generally suggest a bell-shaped distribution over 
the relevant ranges for both ms/subgiant companions and 
horizontal-branch companions, although other distributions 
are not ruled out. To explore a plausible range of pos- 
sibilities, we have studied the following distributions, as 
detailed in Table (a) the above reciprocal distribution 
f{a) oc 1/a, (b) a uniform distribution /(a) — const, (c) a 
linear distribution /(a) oc a, and (d) a Gaussian distribu- 
tion /(a) cx exp[— (a — ao)^/(7^] with appropriately chosen 
central value ao and spread a given in Tabled The ultimate 
purpose is to assess the sensitivity (or lack thereof) of the 
final results on this distribution, as we shall see. 

Calculation of T{"f) clarifies the following points. Pri- 
mordial binaries have a range of a^s whose upper limit is 
considerably larger than that for ms/early-subgiant binaries, 
but most of those binaries which lie between these two up- 
per limits are too wide to be of any practical importance in 
this problem. Thus, it appears that we need consider in de- 
tail only PXBs with ms/early-subgiant companions for our 
purposes here, and Fig. |H] shows the distribution-averaged 
shrinkage time t(7) as a function of the encounter-rate mea- 
sure 7 for such binaries. As can be seen, the above 7"*/^- 
scaling is almost exact at high values of 7, say for 7 > 10*, 
there being a fall-off from this scaling at intermediate 7s, the 
extent of which depends on the case, as shown. We find that 



10' ■ 



f(ai) ~ a, 
((a,) ~ constant 
f(aj) ~ 1/aj 
f(ai) ~ Gaussian 



10"' 



10^ 



10" 



Figure 3. t(7) vs. 7 for PXBs; see text. Curves so normalized 
as to have the same "saturation value" tc = 45 Gyr at low values 
of-y- 



Table 2. Values of 70 obtained by fitting equation illit to 
puted r(7) vs. 7 curves in Fio. 



Type of initial 


Value of 


distribution function 


70 


/(ai) ~ Oi 


8.49 X 10^ 


f{ai) ~ constant 


8.74 X 10^ 


f{ai) ~ 1/ai 


1.21 X 10* 


f{ai) ~ Gaussian 


1.06 X 10* 



the above behavior can be well-represented by the analytic 
approximation 



-(7) 



7*/^ +74/5' 



(12) 



where Aq is a constant which depends on the range 
(a™™,a™"^) (and also on the stellar masses, as explained 
above) , and 70 depends on the above and also on the distri- 
bution f{ai). To illustrate the latter effect, we have given in 
Table |5| the inferred values of 70, where the curves begin to 
deviate from the asymptote, for various distributions in the 
case of tidal capture/exchange binaries. 

We can see the trend that, as the distribution of Ui tends 
to emphasize larger and larger values of at over the permis- 
sible range (as happens in going from a /(fli) ~ to a 
uniform distribution f{ai) ~ const., and further to a lin- 
ear distribution /(a^) ~ ai), 70 decreases. The physical rea- 
son for this is straightforward. Collisional hardening, whose 
rate scales with 7, is dominant at large a (scaling as a^, as 
shown by Eqn. ij^J)- Hence, larger values of ai increase the 
relative contribution of collisional hardening to r, making 
it dominant over a larger range of 7, so that the asymptote 
T ~ Ao'y~'^^^ corresponding to pure collisional hardening is 
followed over a larger range of 7, and so 70 becomes smaller. 
It follows that those distributions which emphasize larger 
values of at will lead to smaller values of 70. 
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Finally, at low values of 7 ( about 10^ - 3 x 10''), the 
following aspect of the low-mass companion's evolutionary 
characteristics enters the picture. The value of r{j) calcu- 
lated in the above manner then exceeds the main-sequence 
lifetime Tc of the companion, a simple, widely-used estimate 
for which is 



13 X 10^ (mc/Me)"^-^ yr. 



(13) 



For a typical low-mass companion with Mc ~ O.GMq, there- 
fore, Tc ~ 45 Gyr. When r(7) calculated as above exceeds 
this value of Tc, what happens is that the companion evolves 
into a giant, and so comes into Roche-lobe contact at a time 
« Tc for essentially all plausible values of a at this point, 
irrespective of the calculated value of r(7). This is formally 
equivalent to the statement that t(7) reaches a saturation 
value of Tc at low values of 7, the change-over occurring at 
7 = 7j, such that r(7c) = Tc- Thus, the computed values can 
be analytically approximated by the prescription that r(7) 
is given by Eq. 1121 for 7 > 7c, and by t(7) = Tc for 7 < 7c. 
This is shown in Fig. 13 where we normalize all the curves 
to the above, common "saturation value" Tc = 45 Gyr. Note 
that the lifetimes of GCs are typically ~ 10 — 14 Gyrs, so 
that, in a given GC, only those PXBs which reach Roche- 
lobe contact within its lifetime would be relevant for our 
purposes. What we have shown in Fig. |3 is the formal be- 
havior of the distribution-averaged r(7) for plausible distri- 
butions of ai. For a given GC, only that range of values of 
which corresponds to Roche-lobe contact within its lifetime 
will go into the specific calculation for it. 



3 EVOLUTION OF COMPACT-STAR 
BINARIES IN GLOBULAR CLUSTERS 

How can we test the above scaling? Since r is not directly 
observable, are there possible signatures that its scaling with 
7 might leave in the observed behavior of the populations 
of compact X-ray binaries in globular clusters? We briefly 
consider this question now and suggest possible answers. 

As remarked earlier, PXBs are formed in GCs primarily 
by tidal capture and exchange interactions in the GC core. 
The rate of the former process is proportional to the en- 
counter rate between single stars in the GC core. The latter 
rate is commonly denoted by F in the literature, and it scales 
as r oc p^r'c/vc with the average core density pc, the veloc- 
ity dispersion Vc of the stars in the core, and the core radius 
Tc. We can describe this as a rate of increase of the number 
npxB of PXBs in the GC which is aiF, where ai is a con- 
stant. In an exchange interaction, one of the members of a 
binary consisting of two normal stars is replaced by a heav- 
ier compact star. The rate of this process is proportional 
to the encounter rate between the above two populations. 
Assuming the population of compact stars in a GC to scale 
with the entire stellar population in the GC, and also the 
population of normal-star binaries in a GC to scale with its 
total population, both of which are normally done, the rate 
of the exchange process also scales with the square of the 
stellar density, and therefore with the above F parameter, 
and we can express it in a similar vein as a2F, where 02 is 
another constant. Thus, we can write the entire formation 
rate phenomenologically as aF, where a = ai + a2. 

After formation, two main processes affect the fate of 



the PXBs. The first is the process of the hardening of 
the PXB to the point of Roche-lobe overflow and conver- 
sion into an XB, as discussed in this paper. This proceeds 
on a timescale tpxb, which means that the PXB popula- 
tion decreases on a timescale tpxb, which we can describe 
phenomenologically by a rate of decrease of npxB which 
is upxb/tpxb- This process may be slightly modifled by 
second-order ones, which are normally ignored. For exam- 
ple, during the above hardening, an exchange encounter of 
the PXB with a single normal GC star heavier than the 
companion in the PXB can replace the latter. 

The second process is the destruction of a PXB by its 
encounter with single stars in the GC core. This can happen 
in the two following ways. First, a fraction of the star-PXB 
encounters leads to a disruption or tonization of the PXB. 
This leads to a reduction in npxB, whose rate is propor- 
tional to 7, the star-binary encounter rate introduced ear- 
lier, and also to npxB- We can thus express this rate of 
decrease in npxB phenomenologically as Pi^fUpxB, where 
/3i is a constant. Secondly, in a smaller fraction of such en- 
counters, a compact star can replace the normal low-mass 
companion in an exchange encounter, resulting in the for- 
mation of a double compact-star binary. This is equivalent 
to destroying the PXB, since such a binary will not evolve 
into an XB. The rate of this process scales with both npxB 
and the rate of encounter between a given PXB and com- 
pact stars. If we again argue that the population of compact 
stars in a GC scales with the entire stellar population in it, 
this rate is ex 7, and the rate of reduction of npxB can be 
written phenomenologically as l32ynpxB, where (52 is a con- 
stant. The total PXB destruction rate can thus be written 
as P'ynpxB , with (3 = Pi + ^2- 



3.1 A Simple "Toy" Evolutionary Scheme 

We now combine the above points into a simple "toy" de- 
scription of PXB and XB evolution in GCs, which we can use 
in an attempt to extract possible signatures of the scaling de- 
scribed in this work. In this "toy" schem e, which is similar in 
spirit to that of lWhite fc GhoshI ^1998^ and lChosh fc White! 
(2001) for following the evolution of X-ray binary popula- 
tions of galaxies outside GCs, the evolution of the PXB pop- 
ulation is given by: 



dnpxB „ „ npxB 
■ al — pjnpxB 



dt 



TpxB 



(14) 



wherein the above rates of increase and decrease of npxB 
have simply been collected together. 

The evolution of the XB population nxB resulting from 
the above PXBs is described in a similar manner: 

dnxB npxB nxB 



dt 



TpXB 



TXB 



(15) 



Here, txb is the evolutionary timescale for XBs. The 
idea here is that XBs are created from PXBs at the rate 
npxB /tpxb , and conclude their mass-transfer phase, and 
so their lifetime as XBs, on a timescale txb ■ 

In the spirit of the "toy" model, all timescales in equa- 
tions 11411 and 11511 can be considered constants, as can a and 
13, while in reality they depend on orbital parameters and 
stellar properties, as also on other parameters. These equa- 
tions can then be solved readily, and of interest to us here 
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is the asymptotic behavior, obtained by setting the time- 
derivatives to zero in these, which yields an XB population: 



nxB = 



(16) 



The effect of collisional hardening is immediately seen on 
the right-hand side of Eqn. llfciL in the second term in the 
denominator. Since collisional hardening always decreases 
T, it increases nxB, other things being equal. This enhance- 
ment in nxB is as expected, as collisional hardening makes 
a larger number PXBs reach Roche-lobe contact. Thus, in 
a GC with given properties, the number of X-ray sources is 
expected to be enhanced by collisional hardening compared 
to what it would be if this effect were negligible. 



3.2 Signature of Collisional Hardening? 

Can we look for observational evidence of the above enhance- 
ment in XB populations of GCs expected from collisional 
hardening? We discuss briefly an attempt to use Chandra 
observations of GCs to this end, with the cautionary re- 
mark that our evolutionary model, as given in the previous 
subsection, is still too simple-minded to apply quantitatively 
to actual GC data. What we are looking for, therefore, is a 
possible qualitative trend that is consistent with the ideas 
of collisional hardening discussed in this paper, which will 
encourage us to perform a more detailed study. 

The trend given by Eqn. 1151 readily translates into one 
of the form 



trend on inclusion of collisional hardening is precisely due 
to the scaling discussed in this paper. However, we stress 
again that ours is only a "toy" model at this stage, relevant 
only for exploring feasibility. To study the effect of ai dis- 
tributions, we have normalized the constants A and B for 
each distribution by having the curve pass through two cho- 
sen points at the lowest and highest values of 7 for which 
data is available. The results show clearly that varying the 
distribution has almost no effect on the trend. 



4 DISCUSSIONS 

The purpose of this paper is to point out an essential effect 
in the hardening of PXBs in GCs, vtz., that collisional hard- 
ening increases with increasing a and orbital period, while 
that due to gravitational radiation has the opposite trend, 
so that their interplay leads to a characteristic scaling of 
the total hardening rate with GC parameters. In our intro- 
ductory treatment of this effect here, we have given very 
simple formulations of many physical processes. First, colli- 
sional hardening is an inherently stochastic process, wherein 
individual events of varying sizes accumulate to yield a final 
state, and the ShuU (1979) rate we have used is a contin- 
uous approximation to it. Secondly, the essential two- and 
three-body interactions that determine the evolution and 
fate of a PXB in a GC are also stochastic by nature. For 
example, approximating an ionization event — in which a 
binary is disrupted — by a continuous term is necessarily a 
great oversimplification. Thus, an improved treatment must 
include a proper formulation of these stochastic processes, 
which we are developing currently. 

Thirdly, we have confined ourselves to circular orbits 
here, while binaries created by tidal capture and/or ex- 
change interactions often have quite eccentric orbits, in 
which tidal circularization must play a dominant role during 
initial phases of hardening. Fourthly, mass segregation is an 
essential effect in GCs, which reflects itself in the accumu- 
lation of the heaviest objects in the core of a GC, and so 
in a change in the effective mass- function in the core. Fi- 
nally, the evolution of the GC must be taken into account in 
a realistic calculation: this would make the GC parameters 
time-dependent, while we have treated them as constants 
here, and may indeed have a large effect if core collapse and 
bounce are modelled. 

The above processes will need to be adequately mod- 
elled in a realistic calculation of the evolution of PXBs and 
XBs in GCs. Such calculations are under way, and the results 
will be reported elsewhere. 

Acknowledgments: It is a pleasure to thank the referee 
for thoughtful comments, which improved the paper consid- 
erably. 



nxB 

where A = I/otxb and B 



(17) 



PA are constants, indepen- 
dent of the cluster parameters. We can compare this with 
d ata obtained f r om C handra observations of GCs, as given 
in lPoolev et"al] ll2003h . This is shown in Fig.|Il 

It is clear that the trend suggested by Eq. I17II is con- 
sistent with the data, while that expected when collisional 
hardening is entirely neglected is not. The flattening of the 
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Defining S = B'y and substituting 5aP = in the above, we 
get 

h^S-ni]l^^, (A3) 
where the indefinite integral I{x) is given by, 
x'^dx 



l+x'i 



(A4) 



Standard expressions for integrals of type / are given in 
Gradshteyn and Ryzhik (1980), from which we get, 

1 



I{x) 



-ln(l + .; ^ 



1 — 2a; cos h 

5 



+ cos — In I 1 + 2a; cos h x 

5 V 5 



(A5) 



. 71" , -1 

sm — tan 

5 



a; - cos ^ \ 2tv _i 

— ' + sm — tan 
5 



From equations lAlL <A2l l. IIA3l l. IIA5II . 

-- -- ~i =b 



(A6) 



APPENDIX A: ANALYTICAL EXPRESSION 

FOR r(7) 

We drop the magnetic braking term in the integral on right- 
hand side of Eqn. (|7J, as explained in the text, and obtain: 

Tpxs(ai,7)« / ^ = ~^^^' (^1) 

Jaf Q^Gwa + aca 7 acw 

where, 

r 3^^. ■ B^J^ (A2) 



